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TRANSCUTANEOUS PETECTION OF ARTERIAL BLEEPING BY TWO- 
PIMENSIONAJ. ULTRASONIC IMAGING QF TISSUE VIBRATIONS 

Qovermnent Rights 

This invention was fimdecj at least in part with grants from the U.S. Office of 
5 N^val Research (ONR) (Award Nos.N00014-96-l-0630, N00014-01-G-0460, and 
NOOO 14-99- 1-0982), and grants from tlie U.S. Department of Defense (DoD) (ARMY 
MRMC Award Nos. DAMD 17-02-2-00 1 4, and DAMP17-00-2-0063), and the U.S. 
government may have certain rights in tliis invention. 

Field of the Invention 

10 This invention generally relates to a mediod for identifying tissue vibrations 

using an ensemble of two-dimensional (2D) ultrasound data acquired with an ultrasound 
transducer, and more specifically, to a method and systerri for identifying bleediag sites 
withiii a body of a patient in real time, based upon tlie tissue vibration that occurs at the 
sites, so that arterial bleeding and bleeding from organs can readily be localized in an 

1 5 image. 

Baclcground of the Invention 

Intemal bleeding is a significant cause of death in cases of trauma, and rapid and 
effective diagnosis of patients with uncontrolled bleeding has long been recognized as an 
important goal to lower mortality and morbidity. Currently, internal bleeding is 

20 diagnosed using angiography. Computed Torjiography (CT), diagnostic peritoneal 
lavage, B-mode ultrasound, and exploratory laparotomy, wliile surgical intervention is 
the most common treatment option. The detection of a bleeding site (i.e., diagnostic 
imaging) and the closure of a bleeding wound should preferably be accomplished as 
quicldy as possible to minimize blood loss in a patient and reduce mortahty and 

25 morbidity associated with such blood loss. 

Advances in duplex and color-flow ultrasound in tlie last two decades have had a 
significant clinical impact on vascular diagnosis. For example, tlie use of Doppler 
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ultrasound has been shown to be effective for targeting a bleeding site, as disclosed in 4 
paper by R.W. Martin, S. Vaezy, P. Kaczkowslci, G. Keilman, S. Carter, M. Caps, K.W. 
Beach, M.I. Plett, and L.A. CxTim, entitled "Hemostasis of punctured vessels using 
Doppler-guided high-intensity iiitrasound," Ulp'asound Med. Biol^ vol. 25, pp. 985-990, 
5 1 999. However, this technique suffers frorn the disadvantage of imaging a limited region 
of interest. Although color-flow ultrasoi^d can image a large region of interest, currently 
it lacks sufficient sensitivity for diagnosing internal bleeding, due to weak scattering from 
blood and the slow flow velocity of blood bleeding from a wound, especially in the case 
of deep bleeds and organ bleeds. In a paper by X. Shi, R.W. Martin, S, Vae2y, and L .A. 

10 Crum, entitled "Quantitative ifl^yestigation of acoustic streaming in blood," J. Acoust 
Soc. Am,y vol. Ill, pp. 1110-1121, 2002, the use of acoustic streaming is proposed for 
distinguishing between stagnant blood and tissue using color-flow UTiages. However this 
paper does not suggest how to detect a bleeding site witli acoustic streaming. The use of 
contrast agents has also shown to be promising for localizing active bleeding sites, as 

15 discussed in a paper by J-B Liu, D.A. Merton, B.B. Goldberg, N.M. Rawool, W.T. Shi, 
and F. Forsberg, entitled "Contrast-enhanced two- and three-dimensional sonography for 
revaluation of intr^-abdominal hemorrhage," J, Ultrasound Med^ vol 21, 
pp. 161-169, 2002. Yet, the use of contrast agents is time-consxmiing and caii sometimes 
be dangerous. 

20 Each of these prior art techniques for detecting bleeding is unable to provide 

efficient real-time images iri which the location of a bleeding site can be rapidly 
identified. Also, as discussed in gieater detail below, srniply using conventional 
color-flow data for imaging a site cannot readily distinguish between pooled blood and 
bleeding at the site. Accordingly, it is necessary to develop a new method for rapid 

25 diagnosis of internal bleeding. 

Physical examination is an important element of tlie initial assessment of a 
trauma patient witli suspected internal bleeding. If an audible "bruit" or a palpable 
"thrill" is found upon physical examination, fiirther diagnostic tests for internal bleeding 
or surgical intervention are often recommended. It has now been established that bruits 

30 and thrills are produced by the forces exerted on yessel walls by eddies created as blood 
flows from a high-pressure region to a low-pressure region through a narrow orifice. The 
pressure fluctuations in eddies cause local vibrations in the vessel wall and surroimding 
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tissue and manifest either as bruits or thrills at the skin surface. The power spectruna of 
the vibrations exhibits a frequency pealc called the "break frequency," which is directly 
related to the diameter of the orifice and tlie local flow velocity through the Sfrouhal 
number. In conventional color-flow ultrasound images, tissue vibrations from abnoniial 
5 bIoo4 flow produce characteristic speckled artifacts in the surrounding tissue. But, these 
artifacts are difficult to interpret and are not quantitative. Tissue vibrations have been 
previously studied using one-dimensional (ID) pulsed Doppler ultrasound, and the prior 
art includes disclosure of a wavelet-based method for detecting and characterizing arterial 
vilDrations in internal bleeding (MI. Plett^ 'TJltrasonic arterial vibrometry with 

10 wavelet-based detectiori and estunation," PhD, dissertation, Univ. of Washington, 2000). 
However, this pulsed Doppler-based technique also has a limited field of view, and is 
along a single scan line. Furthermore, as disclosed in this paper, the processing was done 
off-line, so it was not possible to create images interactively in real time. 

Accordingly, it woiild be desirable to develop a new tissue vibration detectipp 

15 and irnaging rnode for xiltrasound instniments in which vibrations produced by blood 
flow eddies can be detected and color-coded according to their amplitude and frequency 
and overlaid on a B-mode and/or a color-flow image in real time. The tissue vibration 
imaging mode rnight then be used for locating the origin of tlie vibration more precisely, 
relative to the patient's anatomy and/or for obtaining simultaneous information about 

20 vibrations and die underlying blood flow. Acoustic hemostasis using High Intensity 
Focused Ultrasound (HIFU) is a promising new technique for stopping intemal bleeding 
without invasive surgical intervention. For effective targeting and monitoring of 
hemostasis, non-invasive real-tinie localization of a bleeding site in real time is essential. 
This new technique might thus be used both for diagnostic determination of a bleeding 

25 site, and optionally, in conjimction with HIFU or other desired therapy, for localizing the 
bleeding site in real time, so that therapy miglit be effectively applied to stop the bleeding 
as indicated m the U.S. Patent No: 6,083,159, "Metliods and devices for providing 
acoustic hemostasis", U.S. Patent No. 5,993,389, "Devices for providing acoustic 
hemostasis", and US patent 5,882,302, "Methods and devices for providing acoustic 

30 hemostasis". 
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Summary of the Invention 

This invention was developed to localize a bleeding site by imaging tissue 
vibrations caused by arterial bleeding into a free space. Sounds caused by tliese 
vibrations (bruits) are sometimes audible using a stethoscope or palpable at the skin 
5 surface and are indicative of internal bleeding in a traurna patient. Previously, 4 
Doppler-b^ed algorithn^ might have been used for ofiline analysis of vibrations along a 
single xiltrasound scan line. . In contrast, this invention employs an algorithm that 
processes an ensenible of 2D ultrasoiuid data for detecting and imaging tissue vibrations 
in 4 relatively large region of interest. The algorithrn has been implemerited in a 
10 programniable ultr^ound systerp to study the usefiilness of tissue vibrations in real-tirne 
localization of bleeding sites in aiiimals and has been shown to be very effective for this 
purpose. 

The vibration imaging algoritlun uses an ensemble of 2D ultrasoimd data 
acquired during conventional ultiasonic imaging and utilizes the clutter signal, which is 

15 norrnally suppressed in conventional color-flow imaging, to detect and characterize tissue 
vibrations- Thiee signal processing algorithms have been developed that are suitable for 
this pmpose, including one based on phase decomposition, a second based on spectral 
estimation using eigendecomposition, and a tliird based on spectral estimation using 
autoregressive modeling for isolating vibrations from clutter, blood flow, and noise in an 

20 ensemble of 2D ultrasovind data. 

Real-time tissue vibration imaging has been implemented at frame rates, for 
example, of 10 frames/second, ori an ultrasound system with a software-programmable 
signal and image processing back-end. The preliminary results confirm that vibrations 
produced as a result of arterial bleeding can be detected and imaged with this invention. 

25 The vibration amplitude is expected to be the largest near the site of the bleedmg, and this 
fact can be used to localize a bleeding site quicldy and non-invasively. The strong 
backscattered ultiasonic echoes from tissue vibrations can improve visualization of 
interrial bleeding sites that are otherwise hard to image due to weak scattering from 
blood. 

30 Potentially, this new tissue vibration imaging teclinology could be usefiil in a 

variety of devices and clinical settings. For example, a low-cost portable screening 
device with tissue ^db^ation detection flinctionality could be beneficially employed by 
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paramedics and tramna centers to localize bleeding in patients. In addition, a tissye 
vibration imaging mode on high-end ultrasound systems can augment duplex ultrasound 
for enhanced diagnostic capability. The detected tissue vibrations indicative of internal 
bleeding could be presented as an audible signal in a manner recognizable to a person 
5 trained to Usten to bruits using a stethoscope or as a palpable signal recognizable to a 
person trained to detect palpable thrills. 

One aspect of the invention is directed to a method for detecting and localizing 
internal bleeding using an ensemble of 2D ultrasoun4 data by detecting and 
characterizing tissue yibrations caused by blood flow eddies at the bleeding site. The 

10 method includes the step of processing an ensemble of 2D ultrasound data, producing a 
tissue motion spectrum signal at a site. The tissue rnotion spectrum sigiial is then filtered 
to produce a filtered signal from which any contribution to the tissue motion froni a 
source other than bleeding at the site has been substantially minimized, producing a 
tissue vibration signal. A vibration irnage is displayed using the filtered signal and 

1 5 indicates a location of bleeding at the site. 

One approach for processing the ensemble of 2D ultrasound data comprises the 
steps of estimating a correlation matrix from the ultrasound data, and carrying out an 
eigendecomposition of the correlation rnatiix to identify a signal subspace and a noise 
subspace, A frequency of the dominant vibration components in the sigiial subspace and 

20 the noise subspace are then estimated, and based upon that estimate, a vibration 
amplitude estimate and a vibration frequency estimate are determined. At least one of 
the vibration amplitude estimate and the vibration frequency estimate comprises the 
tissue vibration signal. 

A second approach for processing the ensemble of 2D ultrasound data includes 

25 tlie step of computing reflection coefficients from the ultrasoimd data. Linear prediction 
filter coefficients are computed from tlie reflection coefficients. A power specfrum is 
estimated, and the pealcs iii the power spectrum are detected. Based upon the estimate of 
the power spectrum and the pealcs, a vibration amplitude estimate and a vibration 
frequency estimate are detennined. At least one of tlie vibration amplitude estimate and 

30 the vibration frequency estunate again comprises tlie tissue vibration signal. 

In yet another approach for processing the ensemble of 2D ultrasound data, a 
mean clutter velocity is estimated from tlie ultiasound data using autoconelation. Tlie 
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ensemble of 2D ultrasound data are dowr^-mixed with the niean clutter velocity, 
producing a down-mixed sigrial. A phase of the down-rnixed signal and a mean phase of 
the dovyn-mixed signal are determined, md the mean phase is subtracted from the phase 
of the down-mixed signal, producing a residual phase. The residual phase is then 
5 decomposed into its dorninant: cornponents. By applying eriergy and frequency 
thresholds, any contribution to tlie tissue vibration due to noise and blood flow are 
substantially suppressed, yielding an estimate of vibration amplitude and vibration 
frequency of tissue due to bleeding. 

Tlie step of decomposing the residual phase preferably comprises the steps of 

10 estimating a correlation matrix from the residual phase, and performing m 
eigendecomposition of the correlation matrix to determine the dominant components. 

The step of filtering preferably comprises the step of filtering out clutter and 
noise at frequencies that are substantially lower than an expecte4 frequency range of 
tissue vibrations corresponding to bleeding at the site, and also preferably includes the 

15 step of filtering out noise that is at frequencies, which are substantially higher than an 
expected frequency range of tissue vibrations corresponding to bleeding at the site. 

In addition, the method can include the step of confirming that vibrations 
displayed in the vibration image correspond to bleeding at the site, by placing a Doppler 
sample volume at a location of tlie tissue vibration determined from the tissue vibration 

20 image. In this step, the tissue vibration spectriun determined from a relatively larger 
number of samples (typically, 64-512) available using Doppler data can be ernployed to 
confirm the tissue vibration detected from only a few ensemble (typically, 6-16) of 2D 
ultrasoiuid data. 

Tlie step of displaying the vibration image preferably comprises the step of 
25 displaying at least one of a vibration amplitude image and a vibration fiequency image of 
the site. Because of its efficiency, the method can include tlie step of displaying tlie 
vibration image in connection witli an uiiderlying anatomy of tiie site (i.e., the B-mode 
grayscale image), substantially in real time. 

Another aspect of the present invention is directed to apparatus for detecting and 
30 localizing bleeding at an internal site using an ensemble of 2D ultrasound data. Tlie 
apparatus includes an ultrasound transducer for transmitting ultiasound pulses toward the 
internal site and receiving ultrasound data from scatterers at the internal site, including 
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tissue that is vibrating due to bleeding at the site. In one embodiment, a front-end system 
controls tlie iiltrasound pulses produced by the ultrasound transducer and demodulates 
the echoes received by the ultrasound transducer, producing a color-flow signal having 
both in-phase and quadrature components. The apparatus also includes a back-end 
5 system to receive the color-flow signal froni the front-end systern and which includes a 
tissue vibration processor. The fr ont and back-end systems can optionally be combined 
into a smgle unit, or one or piore parts of these systems can be operating remotely from 
other parts of these systems. The tissue vibration processor processes the tiltrasound 
signal to estimate tissue vibrations caused by bleeding, producing a tissue vibration 

10 signal. In one embodiment, the tissue vibration signal is converted to an image signal by 
the back-end system. A display is coupled to the back-end system to receive the image 
signal, to display a tissue vibration image in which bleeding at the internal site is 
indicated. Optionally, the display could be decoupled from the tissue vibration detection 
and identification apparatus. For example, the display can be located physically in a 

15 hospital, while the tissue vibration detection and identification apparatus is located 
physically at another location, such as in an ambulance. Also, instead of a visible 
display, the result of tissue vibration detection and identification can be presented as an 
audible or a palpable output indicating tissue vibrations. The tissue vibration signal can 
also be interpreted by an automated algoritiim to indicate internal bleeding, and tiie result 

20 of the automated interpretation presented as an electronic readout. Generally, the 
fimctions performed by die apparatus aie consistent with the steps of the method 
described above. 

Qrief Pescrlption of the Drawing Figures 

The foregoing aspects and many of flie attendant advantages of this invention will 
25 become more readily appreciated as the same becomes better understood by reference to 
the following detailed description, when taken in conjunction with the accompanyiiig 
drawings, wherein: 

FIGURE 1 is a fiinctional block diagram of an exemplary ultrasoxmd system that 
is suitable for carrying out tissue vibration imaging in accord with the present invention; 
30 FIGURE 2 is a graph of an expected Doppler spectrum (i.e., frequency vs. 

power) from a vibrating sample volume; 
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FIGURE 3 A is a flow chart showing the logical steps of a subspace-based 
algorithm for creating a vibrating tissue image in which bleeding is evident; 

FIGURE 3B is a flow chart shoNving the logical steps of an alternative algorithra 
that uses autoregression for creating a vibrating tissue image in which bleeding is 
5 evident; 

FIGURE 3C is a flow chart showing the logical steps of yet another alternative 
algorithni that uses phase decomposition for creating a vibrating tissue image in which 
bleeding is evident; 

FIGURE 4 is a schematic diagram showing how a gel standoff is used for 
10 acoustically coupling ari ultrasourid transducer to an organ that is exposed and incised to 
create a bleeding site; 

FIGURES 5A, 5B, and 5C are respectively a color-flow image of a punctured 
femoral artery, a vibration amplitude image of the artery, and a vibration frequency 
image of the artery; 

1 5 FIGURES 6A, 6B, and 6C respectively illustrate a color-flow image of an incised 

spleen, a vibration amplitude image of the spleen, and a vibration frequency image of the 
spleen; and 

FIGURES 7A, 7B, and 7G respectively illustrate a color-flow irnage of an incised 
liver, a vibration amplitude image of the liver, and a vibration frequency image of the 
20 organ. 

Descrlplipii Qf the Preferred Embodiment 

Tissue Vibration Imaging System 

FIGURE 1 is a block diagram illustrating an ultiasound system 1 0, which is 
generally similar to a conventional ultrasound system, but which has been modified to 

25 include tissue vibration imaging and is thus usable in practicing the present invention. 
Ultiasound system 10 includes an iiltrasound tiansducer 12 tliat transmits a signal, which 
is modulated with a carrier frequency, typically 1 MHz - 15 MHz, using multiple cycles 
(i.e., 2-20 cycles). The transmitted signal is reflected by scatterers (not shown) along the 
beain path and is received after a time delay, which depends upon the distance of the 

30 scatter from the transducer. In an acquisition stage, the acoustic echoes received fr om the 
tissue are converted to electrical sigiials by the transducer, and these signals are digitized 
by analog-to-digital converters (not separately shown), A front-end subsystem 14 
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includes a beam former 16 that performs dynamic focusing, apodization, and steering of 
both transmitted and received ultrasonic waveforms. Also included in front-end 
system 14 are a tinie- gain-compensation (TGC) circuit 18 that amplifies signals vvith a 
variable gain that is proportional to the depth within tissue, and a radio frequency (RF) 

5 demodulator and decimator20 t^iat digitally removes the high frequency carrier by 
quadrature demodulation and decimation, providing both in-phase (I) and quadrature (O) 
samples, which may be represented as a complex quantity, I(t) + jQ(t), The acquired 
quadrature ensemble (or color-flow) data are then processed in a back-end subsystem 22, 
depending on the one (or more) ultrasoimd mode(s) that is/are selected, e.g., B-mode, 

1 0 color-flow mode, tissue vibration mode, and Doppler mode. 

For producing anatomic jinages of tissue, the signal of interest is the envelope of 
I(t) + jQ(0' A B-mode processor 24 computes tlie magnitude of the echo, 
5a (0 = (0+2^(0 and compresses the dynamic range to make it suitable for display as 
a grayscale image on a monitor. The time delay introduced by the scatterers is reflected 

15 in the phase of the complex quantity I(t) H- jQ(t), Thus, the phase of the complex 
received signal provides an estiinate of tlie instantaneous position of the scatterer. By 
monitoring the change of phase over tirne, the displacement and velocity of the scatterer 
can be estimated. In color-flow imaging, multiple pulses (commonly from 6 to 
16 pulses) are transmitted and received along each scan line at a rate known as the pulse 

20 repetition frequency (PRF). A collection of received temporal samples from each spatial 
location is tluis called an "ensemble." A color-flow processor 26 estimates the blood 
flow velocity from the ensemble of data by estimating the phase difference between the 
adjacent temporal samples, typically using an autocorrelation algorithm. A 2D inaage is 
created by acquiring multiple samples from different spatial locations. In Doppler mode 

25 which is implemented by a Doppler processor 30, scanning is performed along a single 
scan line, and a spectrum of the blood velocity from a single spatial location is estimated 
from a substantially larger ensemble of data (typically, data from 64-512 pulses). Before 
displaying the processed image frame on a raster monitor or display 36, scan conversion 
is performed by a scan converter circuit 32, which converts the acqiured ultrasound data 

30 from polar coordinates to the Cartesian coordinates used by the raster display. Post 
processing may optionally be appUed by a post-processing circuit 34, to improve the 
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quality of the displayed image, as well eis to combine ttie anatomy and flow images on 
the display. 

A tissue vibration processor 28 that is used to process tlie ultrasound data in the 
present invention is shown iix FIGURE 1. The ultrasound data acquisition for tissue 
5 vibration processing and imaging is tJre same as that used for color-flow, so that 
quadrature data eriserable or color-flow data aie input to the tissue vibration processor. 
Ho weyer,. instead of estimating blood flow velocity from these data, the tissue vibration 
processor estimates the instantaneous displacement of the scattering tissue fi:om the phase 
of the complex received signal. This tissue motion is referred to as clutter in 

10 conventional color-flovv iinaging and is suppressed using clutter filters. Typically, 
cardiac pulsation, respiratioji and transducer motion each cari contribute to an observed 
displacement or motion of tissue. Si}ch motion is at a low frequency of a few Hertz or 
less. When iiitemal bleeding is present, the tissue vibrates locally vyith a frequency 
ranging from a few tens of Heitz to a few hundred Hertz. By analyzing the frequencies 

15 of tlie different cornponents of tissue motion, vibrations caused by bleeding may be 
distinguished from clutter caused by other sources of movement The tissue vibration 
processor performs this analysis by decomposing the tissue motion into the dominant 
motion components and identifying any motion components that appear to be at a 
frequency higher than that of caidiac pulsation. 

20 It is contemplated that tissue vibration processor 28 can be implemented as an 

additional fixed-ftinction circuit board or application specific integrated circuit (ASIC) in 
conventional ultrasomd machines. Optionally, the tissue vibration processor can be 
combined with color-flow processor 26, since both process tiie same data ensemble. A 
standalone tissue vibration unaging device can be implemented with front-end 

25 subsystem 14, B-mode processor 24, tissue vibration processor 28, and scan 
converter 32. Those of ordinary skill in the art will appreciate that the tissue vibration 
processor can be implemented in software/hardware using one or more digital signal 
processors (DSPs) or alternatively, in an ASIC, or even on a conventional general 
purpose processor chip that access niachuie language instructions stored in a memory 

30 accessed by the processor to c^^ry out tlie processing steps of the tissue vibration 
processor. 
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The computational power of ultrasound machines has increased significantly in 
recetit years, benefiting fi:om advances in processor technology. Thus, the additional 
computational burden of the tissne vibration imaging algorithms discussed below can be 
reasonably supported in modem ultrasound machines. Previously, a programmable 
5 ultrasound signal and image processing system suitable for vise as the tissue vibration 
processor was developed using a new generation of high-performance multimedia 
processors to support all of the conventional processirig modes, such as B, M, color-flow, 
and Doppler in software (Sikdar S, Shamdasani V, Gong L, Managuli R, Hayashi T, 
Mitake T, Kim Y. "A single mediaprocessor-based programTnable tiltrasound system" 
10 DSEE Trans Inf. Tech. Biomed 2003; 7:64-70), and subsequently, was shown to be useful 
in implementing tissue vibration processing in accord with the present invention. The 
main strength of a programmable system is the ease of development of new modes and 
applications such as tissue yibration unaging without the need for hardware 
modifications to conventional ultrasound machines. Integrated tissue vibration imaging 
15 using the software-programmable ultrasound system has tlius been effectively and 
beneficially used for realrtime visualization of vibrations in 2D ultrasound scans. 
Algorithms for Tissue Vibration Imagine 

In conventional color-flow imaging, tlie velocity of blood flow is estimated by 
computing the average phase difference between multiple ultrasound echoes (typically 
20 6-1 6 pulses) tliat are received firom a sample volume. Echoes firom moving tissue tend to 
have a significantly higher signal strength (typically 40 dB-60 dB higher), compared to 
the wealc scattering from blood, and also have lower velocities. This high amplitude and 
low frequency tissue signal is commonly referred to as clutter aiid tends to bias the 
estimated blood flow velocity. Thus, clutter is suppressed using appropriate filters in 
25 conventional color flow imaging. The main components of clutter are cardiac pulsation, 
respiration, and transducer movement When blood flow eddies are present, any local 
tissue vibrations, e.g., those caused by the blood flow eddies at a bleeding site, will also 
be part of this clutter and would normally be suppressed in conventional ultrasound 
processing systems. 

30 Tlie present invention separates the tissue vibrations from the remaining clutter 

and flow signals. In achieving tliis function, it was recognized tliat the tissue vibrations 
and clutter produce statistically independent signals tliat have different frequency content. 



\VO 2005/037060'*' 



-12- 



PGT/yS2004/032427 



While clutter due to cardiac pulsation and breathing typically occurs at 1 Hz or less, 
tissue vibrations typically occur at 50 Hz or more. Other noise sources are at 
substantially higher frequencies. Scattering from tissue is typically more coherent 
compared to the scattering from blood, because the tissue scatterers are more tightly 

5 bound together and tend to move as a group. Thus, compared to the clutter from other 
sources ^nd tissue vibration signals, tlie blood flow signal typically has a much greater 
frequency bandwidth. Due to its weak signal strength and greater bandwidth, blood flow 
signals may be considered as noise compared to the stronger and rnore coherent tissue 
vibration sigrials for purposes of tliis invention. Therefore, tissue vibrations can be 

10 distinguished from clutter an4 blood flow based on spectral analysis. Spectral analysis of 
the phase of the received ultrasound echo can be used to separate the components of tlie 
scatterer motion, ignoring the scattered signal strength, whereas spectral analysis of the 
complex ultrasound echo considers both the signal strengtli and the motion components. 

Due to the limited number of temporal ultrasound samples (6-16 pulses) 

15 preferably used in the present invention, conventional Fourier transform and filtering 
approaches lack sufficient resolution for separating the signals. Tlrerefore, this invention 
employs high-resolution specfral estimation techniques to carry out this fimction. Two 
high-resolution spectral estimation techniques were identified as suitable for this purpose, 
including eigendecomposition-based spectial estiination, wliich models the signal as an 

20 optimum set of orthogonal components, and autoregressive spectral estimation, which 
models the signal as the output of an autoregressive linear prediction filter driven by 
white Gaussian noise. Accordingly, three signal processing algorithms were developed 
for isolating tissue vibrations. The first algorithm is based on an 
eigendecomposition-based spectral analysis of the phase of the receiyed ultrasound echo; 

25 the second algorithm is based on an eigendecomposition-based spectral analysis of the 
complex ultrasound echo; and, the tliird algorithm is based on an autoregressive specfral 
analysis of tlie complex ultrasound echo. Since eigendecomposition is a 
computationally-intensive operation, an approximate eigendecomposition u tilizin g 
iterative QR factorization is use4 to develop a computationally-efficient algorithm. 

30 Signal Model 

To model the received signal firom vibrating tissue, the tissue being imaged is 
approximated with S point scatterers having uniform motion and randomly distributed at 
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locations (r,=[r,]e,+[v/,.]e,+[«>..]^,), . = 1-5, in a sample volume where 
denote tlie unit direction vectors in spherical coordinates. The instantaneous 
positioii of the scatterers, v(F, /) , is given by: 

5 where CO. (/). (/)) denote the displacement as a function of time. If the scattering 
from the sample volume is uniform with a as the average scattering coefficient, then the 
scattering fimction of the sample volume isov(r./) • The complex received signal from 
the sample vohiine, j;(t, /) > can then be modeled as a convolution of the pulse echo 
spatial iinpulse response, /i^.(F.r) . of a single point scatterer, the temporal response of 

10 the transducer, and the scattering fiinction, ccv(r,t) . 

yir, 0 = (r, t) * x(t) * av(r,t) + t) 

where the time indices rand t refer to "fast" time and "slow" time, respectively, /o is the 
center frequency of the traiisducer, and niT,t)is white thermal noise. Combining 
Eqs. (1), (2). and (3) results in: 

where c is the speed of sound and ^j^^ two-way pulse propagation time 

between tlie transducer and each point scatterer. The scatterer displacement for 
vibrations is smaU compared to the spatial size of the pulse echo spatial impulse response 
and the envelope of tlie transducer response. Thus, tlie "slow" time variations in the first 
20 two terms niay be neglected and Eq. (3) can be simplified to: 

= A{T)e ' +n(T,t) 

where A(r) is tlie complex amplitude of the scattered signal. It is apparent that the 
complex received signal is phase modulated witli die instantaneous radial displacement. 
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If the scatterers in the sample volume are all vibrating radially in a siinple 
harmonic fashion with peak displacement ao and frequency fyib, the tissue displacernent 
due to cardiac pulsation, breathirig, and other tissue movemerit (i.e., clutter or noise) 
relative to the transducer is dtissit)- This motion will hereinafter be referred to as the 
5 "clutter rnotion " Then, the corribiiied displacement can be considered to be a 
superposition, as follows: 

dr (0 = ^im. (0 + ^0 sin(27/,,,0 (5) 

An ensemble of ultrasound pulses is transmitted in the same direction at a rate 
known as the pulse repetition frequency (PRF). Then, the complex received signal from 
10 the 7«^'' pulse transmission, j<T,m), is: 

y(T,m) = Air)e " ^ +«(r,7;7) 

where Tprf is the pulse repetitiori interval. The Fourier transform of the phase-modulated 
complex received signal is a Bessel series: 



Y(T,f) = A(r)i!3ie - ) 



= ^(r)|c(/) * 'gv V,. ifi)3if - 2m,T,^f^.„ )| + Nir, /) 

/|s-ao 

1 5 where J/ are Bessel ftmctions of the first Icind, ji = — , 5 is the Dirac deUa function, 

c 

and c(J) is the spectrum of the clutter motion (the clutter spectrum), and N{r,J) is the 
noise spectrum. 

FIGURE 2 illustrates a typical power spectrum 40 of the ultrasound signal when 
a tissue vibration is present. The spectrurn includes multiple copies of the clutter 
20 spectrum separated by the vibration frequency, as indicated by Eq. (8), wliich is 
presented below. A low frequency peak 42 at fti^s corresponds to die clutter spectrum, 
while symmetric pealcs /, and /*- indicated respectively by reference numbers 44 and 46 
conespond to vibration, and pu^sy p+y and p., aie tire conesponding peak powers. Tlie 
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frequency peaks at /t and / are referred to herein as a "matching pair." For 
small-amplitude vibrations, higher-order terms can be ignored; thus, most of the spectra] 
energy will be present in the three frequency peaks, /, and respectively. Since, 

vibration ampUtude. Therefore, the vibration frequency and ampUtude may be estimated 
from the power spectrum as follows: 



» /? , the ratio of the power in the frequency peaks can provide an estimate of the 



power 



/.-/- 




(8) 



These estimators are referred to herein as the ''specti'al frequency estimatof^ and 
the "poTver ratio amplitude estimator.^'' Alternatively, tlie vibration frequency and 
10 amplitude may be estimated from tlie residual phase {j^>(A:)}^^^of the ultrasound signal 

after suppressing the effects of clutter motion. A coarse computationally-efScieiit 
estimate of the frequency of the dominant components, /^.^^ , can be obtained by counting 

the zero crossings, iVze/o, ^ the residual phase. This estimate can be fiirther refined by 
interpolating the residual phase to compute the mean period of oscillation The vibration 
15 amplitude may be estimated by the variance of the residual phase. These estimators are 
defined as follows: 



T phase 

Jvib "■ 



/.-/- 



C 



52"^'=7tr^^(<^(i^)) (9) 



and are respectively referred to herein as the ^^zero-crossing frequency estimator*^ and the 
^^phase variance amplitude estimator. 

20 For real-time tissue vibration imaging, only a shoit ensemble of ultrasound data 

(typically, 6-16 pulses or echoes) from each sample volume in a region of interest is 
available for processing. Conventional color-flow imaging systems utilize clutter 
filtering to suppress the clutter, while retaining the blood flow. However, due to the 
small number of temporal samples, tlie conventional clutter filtering-based approach, o^ a 

25 Fourier-based approach laclcs sufficient resolution to discriminate between the tissue 
vibrations, blood flow, and clutter. A paiametric approach that utilizes the characteristics 
of tlie vibration signal appears better suited to make tiiis distinction. Three parametric 
approaches may be taken, based on tlie model of the ultrasound signal in Eqs. (6) and (7), 
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including: (a) estimation of a pair of complex exponentials in noise; (b) autoregressive 
modeliiig; and, (c) decomposition of the phase of the ultrasound signal. In the following 
section, vibration detection algorithm based on these three parametiic approaches are 
described in greater detail. One method of producmg the complex ultrasound signal in 
5 Eq. (4) is a quadrature demodulation of the received ultrasound signal. An alternative 
method is to . compute the time delays producing die phase variations in Eq. (4) by 
processing the received RP ultrasound data using a cross correlation technique. 
Vibration Imaging Using Estimation of Complex Exponentials in Noise 

Using the inverse Fomief transform of the Bessel expansion in Eq. (7), Eq. (6) is 
1 0 expanded, as follows: 

Li>«o J 

Tlius, the ultrasound signal can be modeled as a sum of complex exponentials 
embedded iri noise. As can be seen from the expression enclosed by square brackets In 
Eq. (10), vibrations correspond to matching pairs of complex exponentials. In contrast, 

15 the complex exponentials coiresponding to clutter motion will typically not have such 
matcliing pairs of frequencies. The frequencies (ftiss>f+>f) may be estimated using die 
root-MUSIC and ESPRIT algoriflxms (disclosed by P. Stoica and R. Moses in 
''Introduction to Spectral Analysis,'' Upper Saddle River, NJ: Prentice-Hall, 1997). 
Vibrations may tlien be detected using a matching peak criterion 

20 |/^ + /. - 2/,^,,. I < , arid tiie vibration amplitude and frequency can be estimated 

using Eq (9), which is set forth above. The steps of the algoritiim are described in more 
detail below, in regard to FIGURE 3 A. Based on this criterion, vibrations can be 
detected and distinguished from clutter motion. Any blood flow signals may be 
considered as part of the noise spectrum. 

25 FIGURE 3 A illustrates a flow chart 50 that shows tiie logical steps involved in a 

first algoritiim for estirnating tiie tissue vibrations based upon a pair of complex 
exponentials in clutter or noise, that are normally excluded from color-flow processing. 
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The procedure begins with a quadrature-deiiciodulated ensemble of 2D ultrasound 
data 52. 

In a step 54, a correlatio|i matrix is estimated from the color-flow data. In a 
step 56, the conelation matrix is employe4 to cany out an eigendecompositioq, 
5 producing a signal subspace 58 and a noise subspace 60. Using the signal subspace and 
the noise subspace, the frequency of the dominarit components is estimated in a step 62, 
by employing the root-MUSIC and ESPRIT algorithms, as noted above. A decision 
step 64 then determines if each dominant component is double-sided, while a step 66 
estimates the vibration amplitude and frequency of each dominant component using 

1 0 Eq. (8). If a dominant component is not double-sided (i.e., is not a potential tissue 
vibration component), decision step 64 returns a "zero," while if the dominant 
component is double-sided, the decision step return a one. A multiplier 68 then 
rnultiplies tlie output of decision step 64 times the vibration amplitude and frequency 
estimate for tlie dominant component, yielding a null if the dominant component is not a 

15 tissue vibration component and returning tlie estimate of vibration amplitude and 
frequency of the dominant component otlieiwise. A median filter 70 then filters isolated 
falsely-detected vibrations and other undesired noise from the results, so that the 
remaining vibration image indicating a bleeding site is displayed in a step 72. 
Vibration Imaging Using an Autoregressive Simal Model 

20 The ultrasound signal from vibrations can be modeled as tlie ou^ut of a p^-ovdGv 

autoregressive linear predictioii filter with white Gaussian noise having a variance o^, as 
the input, as follows: 

yir, m) = ^ a,„_^ ir)y(T, m - fc) -V- n(T, m) 

(11) 

Tlie hnear prediction coefficients, ak(r)^ can be computed using either a 
25 least-squares miiumization of the prediction enors or using the computationally-efficient 
Burg algorithm, as explained by Stoica and Moses in the above-referenced paper. A 
higli-resolution spectrgd estimate can then be obtained from this autoregressive model as 
follows: 
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1 + 



J2nkf 



(12) 



From Eq. (8), the presence of symmetric matching pairs of frequency peaks in the 
power spectrum around the clutter motion peak may be detected as a vibration. As 
before, any flow signals may be regarded as noise. For ultrasound signals from 

5 vibrations, die power spectrum in Eq. (13) will have frequency peaks (ftiss»f+yf) at the 

p 

local minima of the polynomial A{r, /) = 1 ^ ^a- (r)g"^^'^ . Vibrations can be 

detected using a matching peak criterion |/^ + /_ - 2 f^^^ | < J^f,,resfioid > and the vibration 

ampUtude and frequency can be estimated using Eq. (9). The steps of this algorithm are 
described in more detail below, in conriection with FIGURE 3B. 

10 As shown in a flow chait 80 in FIGURE 3B, the second alternative algorithm 

also begins with quadrature-demodulated ensemble data set 52. In a step 82, reflection 
coefficients are computed for e^ch ensemble of the quadrature-demodulated data. Using 
the reflection coefficients, linear prediction coefficients are determined in a step 84. In 4 
step 86, the power spectimn is estimated from the liaear predictior^ coefficients and die 

1 5 peaks in the power spectaim are detected. Again, decision step 64 determines if tlie 
pealcs thus identified are for tissue vibration by determining if they are double-sided and 
returning a zero if not, and a one, if so. Also, step 66 provides for estimating the 
Adbration amplitude and frequency at each of these peaks, and the results from decision 
step 64 are multipUed the estimated amplitude and frequency iri multipUer 68. Median 

20 filter 70 is then appUed to the results and tlie filtered image data are displayed as a 
vibration image, in step 72. 

Vibration Imaging Based on Phase Decomposition 

Another model of vibrations can be based on the phase ulti-asound signal. As 
shown in Eq. (7), vibrations vsdll produce an oscillatory signature in the phase, which will 
25 typically not be present in the case of clutter motion, Altliough flow signals may have 
oscillatory phase, the echoes fiom vibrating tissue aie expected to be more coherent tlian 
those fiom flow. Thus, their phase may be modeled by a smaller number of dominant 
components. Accordingly, a vibration detection algorithm can also be based on 
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decomposition of the phase of the ultrasound signal into its dominant components and 
testing for oscillatory phase. Altema.tively, instead of using quadrature-demodulated 
ultrasound data, the phase cari be estimated frona RF yltrasound data by estiinating the 
time delays between a pair of RF ultrasound data.. 
5 Any linear tinie-varying motion is first suppressed by down mixing the ensemble 

of 2P ultrasoiind data with the mean clutter velocity, estirnated using the corxventional 
autocorrelation method. Tlie phase of the ensemble of 2D ultrasound data is then 
cornputed, and the mean phase is subtracted to suppress the effect of the stationary echo. 
The residual phase is then decomposed into its dominant components using ^ metliod 

10 similar to principal cornponent analysis. The first step of the decomposition involves the 
estimation of the con-elation piatrix of the residual phase using the modified covariance 
method (Marple, 1987). An approximate eigendecomposition can then be performed 
using iterative QR factorization of the correlation matrix. Tlie approximate eigenvalues, 
rpay be estuxiated by the diagonal elements of tlie upper triangular matrix Rk after the 

15 A**^ iteration. The eigenvectors are arranged ixi order of decreasing eigenvalues. The 
eigenvalues are a measure of the signal energy contributed by the corresponding 

eigenvector. Tlius, the fi-action of tlie total signal energy contained in the p dominant 

p 

components can be estimated using = — . Therefore, noise and blood flow can 

be suppressed using a threshold criterion, Ep>E threshold' To fiirther separate tissue 
20 vibrations fi-om clutter motion, the fact that tissue vibrations have a higher frequency 
compared to clutter motion is applied- Vibrations can then be separated from clutter 
using a frequency threshold criterion > F,f,resf,ou > where Fthreshoid is chosen so tliat at 
least one half of one period of tlie vibration is contauied in an ensemble. Tlie vibration 
frequency and amplitude may be estimated using Eq. (10). The steps of this algorithm 
25 are described in more detail below, in connection with FIGURE 3C. 

A flow chart 90 illustrates the logical steps of the third algorithm. Again, starting 
with the quadrature ensemble data, a step 92 down mixes the quadrature data with a 
mean clutter velocity, which is determined from the data. A step 94 provides for 
computing an unwrapped phase of the quadrature ensemble or color flow data, and then 
30 subtracting tlie mean clutter velocity from the unwiapped phase, resulting in a residual 
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phase. As explained above, instead of determining the phase froni quadrature- 
demodulated data, the phase can be deterrnined from RF ultrasound data by estiinatipg 
time delays between a pair of RF iiltrasound data. Using the residual phase, a step 96 
estirnates a correlation rriatrix, which is tlien used to carry out a QR factorization ii; a 
5 step 98, yielding an eigenvalue estimate 100, an eigenvector estimate 102, and a vibration 
ampUtifde and frequency estimate 104, which are determined using Eq. (10), as noted 
above. Using the eigenvalue estimate, a decision step 106 determines if the total energy 
contained in tlie p dominant cornponent is greater that a predefined threshold, T. If so, 
decision step 106 returns a zero if not, and a one if so. Similarly, a decision step 108 

10 determines if the estiiiiate eigenvector has a frequency that is greater than a predefiiied 
threshold, F. If so, decision step 108 returns a one, and if not, a zero. The results of 
decision steps 106 and 108, and the estimated vibration amplitude and frequency of the 
dominant components are then multiplied together by a multiplier 1 10, so that if either of 
the decision bloclcs has returned a zero, the result is null, but if neither has retumed a 

15 zero, the estimated vibration amplitude and frequency from step 104 are retumed. Again, 
median filter 70 is applied to tlie estimated amplitude and frequency, providing filtered 
results that are displayed as the vibration image, indicating a bleeding site, in step 72. 
Detectable Vibration Amplitudes and Frequencies 

In experiments using a physical phantona model, tissue vibrations with a peak 

20 amplitude of about 1 |am have been accurately detected. Tlie minimum detectable 
vibration amplitude depends upon die noise level and dynamic range of the phase of the 
receiyed ultrasound echo. In modern ultrasound machines, the phase caii have a dynaniic 
range of 96 dB or more (for 16-bit quadiature-demodulated data) and tlie signal exceeds 
the electronic and thermal noise level by typically 80 dB or more. Tlierefore, frojn 

25 Eq. (4), vibrations as sm^ll as 50 nm may theoretically be detected using ultrasound. 
Practically, tiie attenuation of the ultrasound signal will reduce the dynamic range and 
limit Hie minimum detectable arnplitude in deep tissue to -0.5 \.im. 

The detectable vibration fiequencies depend upon the choice of PRF, i.e, on Fpjar, 
A PRF that is too low compared to the vibration frequency would lead to aliasing, while 

30 selecting a PRF that is too high will fail to detect low-frequency vibrations. A vibration 
can be detected only if at least half of one vibration cycle is captured within the temporal 
window coiTesponding to an ensemble. Tlius, all vibrations with frequency between 
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jl>5£L and —^^can be detected theoretically without aliasing for an ejisemble size E. 

I'^E 2 

Since vibrations can be broadband, a high-frequency vibration interrogated at a low PRF 
value can be mistaken for noise ming this algorithrn. Thus, for better sensitivity, it is 

desirable to select a PRF and an ensernble size so that only a few periods of the vibration 

kF 

5 are included in an ensemble. Accordingly, tlie maximum detectable frequency is — 

E 

when k periods of tlie vibration are included in an ensemble. A simulation and phantota 
experiments that were carried out indicate that reliable detection may be performed using 
only one half to six vibration periods during the interrogation period. For example, with 
a PRF of 1 l<Hz and an ensemble size of 16, vibrations with frequency between 3L3 Hz 

1 0 and 375 Hz may be reliably detected. 
Determiaation of Bleediag Rate 

Since the tissue vibrations are produced by the blood flow eddies, the frequericy 
of the tissue vibrations is the san^e as tlie frequency of the eddies. Hie Strouhal number 
(S) relates the frequency of eddies (fyit) produced by an orifice to the diameter of the 

1 5 orifice (D), and the flow rate through the orifice (U), according to : 

S^IxiL^ (13) 
U 

Since tlie Stioulial number remains constant, tlie vibration frequency can be related to tlie 
diameter of the orifice and tlie flow rate througli the orifice. In addition, the energy in the 
eddies (E) and thus, the amplitude of the tissue vibrations (a^itX is directly proportional to 
20 the flow rate, as follows: 

Eocal,<x:U^ (14) 

Tlie bleeding rate Q = can thus be detemiined fiom the frequency of the 

vibrations and tlie amplitude of tlie tissue vibrations. 
Sources of Artifacts 

25 In color-flow data acquisition, interrogation along each scan line is performed for 

only a brief period of time. Vibrations are transient, with typical durations of 
10 ms-100 ms. Thus, tliere is a possibility that some vibrations may not be interrogated. 
Since tlie vibrations typically have a relatively large spatial extent and repeat every 
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cardiac cycle, it is unlikely that the vibrations will be missed entirely; however, the 
spatial extent of the vibrations visible in the image may be only a part of the true spatial 
extent. By appropriately choosing the PRF and the region of interest, such discrepancies 
may he minimized. 

5 Other artifacts may be falsely detected as vibrations. Transducer motion rnay 

introduce additional frequency peaks in the clutter spectrum and may cause false 
detections; however, these false detections may be minimized by using a trained 
sonographer to perform the scanning. Vibrations in tlie tensed skeletal muscle of the 
sonographer, and any ambient vibrations may be detected in the vibration image. In 

10 addition, the liigh-resolution spectral estimation methods may produce spurious peaks 
that can be falsely detected as vibrations. Such artifacts can be easily distinguished from 
pathological vibrations, which ore expected to be correlated with the anatomy and 
periodic with every cardiac cycle. These artifacts can be also avoided if additional 
temporal samples are available. Any vibrations displayed in tlie vibration image should 

1 5 therefore be confirmed with the vibration spectrum by placing a Doppler sample volume 
at the location of the peak intensity. 
Comparison of tlie Algoritlims Derived fr om Modeling 

Tlie ability of the proposed algorithms to detect vibrations was evaluated using a 
simiJation model. Simulations show that subspace-based algoritlmas such as MUSIC 

20 and ESPRIT have high sensitiyity (96%) and specificity (98%) for detecting narrowband 
vibrations in tlie presence of clutter as well as blood flow and are robust even when 
broadband vibrations are present. For narrowband vibrations, an algorithm based on an 
autoregressive model has a slightly improved specificity (99%), a comparable sensitivity, 
and is robust to broadband vibrations. The phase deconiposition-based algorithm has a 

25 slightly lower sensitivity (93%) and specificity (98%), but is more robust to broadband 
vibrations. 

The computational requirements of tlie proposed algorithms are shown below in 
Table 2. Tlie subspace-based algorithms (MUSIC/ESPRIT) have a computational 
requirement that is highly dependent on the choice of model order. In these algorithms, 
30 eigendecomposition is tlie most computationally-intensive task. The 
autoregression-based algorithm is less computationally intensive and the computational 
requu ement is less dependent on the model order. In this case, tlie computation of the 
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FFT for spectral estimation is tlie most computationally-intensive task. The phase 
decomposition method is the least computationally intensive, since it involves operations 
on real signals only. Tlius, the phase-4ecomposition algorithm is most suitable for 
real-time implementation. 

Tibbie 1 - Field 11 Siinulation Parameters 



Center frequency 
PRF 

Transducer excitation 
Transducer impulse response 
/ nimiber 

Number of elements 
Transducer heiglat 
Element pitch 
Element kerf 
Transmit aperture 
Receive aperture 
Transmit focus 
Receive focus 
Elevation focus 
Matlieraatic^l element size 
Samplmg frequency 
Sound velocity 
Number of scan lines 
Number of ensembles 



5 MHz 

500 Horizontal 
5-period sinusoid 

Hanning- weighted 2-period sinusoid 
2 

192 

15 nun 
0.4 mm 
0.03 mm 
25.6 ram 
25.6 mm 
40 mm 

30 mm to 100 mm in steps of 10 mm 

20 mm 

0.37mm x 1.5 mm 

105 1V1HZ 

1540 m/s 

32 

10 



Table 2 ~ Computational requirement (million operations/sec) for real-time 



imagmg at 10 frames/s vyith 32 scan lines, 256 samples/scan line, and ensemble 10 



Algorithm 


Model Order I 


p=2 


P=3 


p=4 


MUSIC 




3631 


7653 


ESPRIT 




2218 


6455 


AR 




1606 


1630 


Phase-decomposition 


181 


489 


1107 


Color flow 


89 
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In Vivo Empirical Tests of Arterial and Organ Bleeding in Animal Model 

FIGURE 5A illiastrates a color-flow image of a punctured femoral artery. A 
dotted line box indicates the region of interest. An arrow indicates the approxixc^ate 
location of a puncture. In FIGURE 5B, a vibration amplitude inaage of the punctured 
5 femoral artery is illustrated. Although not visible m color in this grayscale image, the 
vibrations at the puncture site in the original color image appear green relative to the 
shades of gray in the surroynding tissue. Similarly, in FIGURE 5C, a vibration 
frequency image shows low-frequency vibrations in color at the puncture site, relative to 
shades of gray for the surrounding tissue. 
10 FIGURE 4 is a schematic view 120 that illustrates how ultrasound transducer 12 

is ernployed to produce color-flow dat^ for an incised organ 124. To couple acoiistic 
signals from and to the ultrasound transducer, a gel standoff 122 is disposed between tlie 
ultrasound transducer and the exposed organ. The gel standoff provides an efficient 
couple between the tissue of the organ and the transducer. In this manner, tissue 
15 vibratior^s caused by bleeding tlirpugh m incision 126 can readily be produced in accord 
witli the present invention. 

In FIGURE 6A, a color-flow image of an incised spleen is shovm. An arrow 
indicates die pooled blood from an incision, where the blood pool has accumulated 
between the organ and the gel stand-off. Sorne blood flow is observed at the bleeding 
20 site as indicated by an arrow. FIGURE 6B shows tlie vibration amplitude image, where 
vibrations surrounding tlie bleeding site are indicated by an arrow. A vibmtion frequency 
image in FIGURE 6C shows die low-frequency vibrations in color to indicate the 
bleeding site. 

FIGURE 7A shows the color-flow image of an incised liver in a rabbit The liver 
25 was surgically exposed and a triangular mcision was made with a scalpel, causing 
moderate bleeding. Tlie liver was imaged tliough a gel standoff. The location of the 
incision is clearly seen in the B-mode image (not shown), and some flow is observed m 
the artery feeding the bleed. However, from tlie color-flow image, it is not possible to 
identify the bleeding site at the incision. Tlie vibration amplitude image is shown in 
30 FIGURE 7B. Vibrations aie clearly observed surrounding the incision from which 
bleeding is occurring, as indicated by an aiTOW. Tlie vibration frequency image in 
FIGURE 7C shows low-frequency vibrations, and the bleeding site is again clearly 
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yisible (although not so apparent in this grayscale image as in the original color images). 
The vibration image and color-flow images therefore provide complementary 
information and a clear correlation between the spatial location of the vibration, and tlie 
underlying anatomy can be observed. Thus, a combination of B-mode, color-flow, and 
5 tissue vibration images can be a powerful diagnostic tool for detecting the location of 
internal bleeding. 
Summary 

Empirical evidence demonstrates tlie feasibility of real-time ultrasoxmd imaging 
of low-intensity local vibrations in the vessel wall and smrounding tissue associated witti 

10 punctured arteries and organ bleeds. Several algorithms based on parametric signal 
decomposition and spectral esttrnation haye been developed for imaging small-amplitude 
tissue vibrations using as few as 10 temporal samples. The vibration amplitude and 
frequency can be estimated accurately, and real-time tissue vibration imaging has been 
impleinented on an xiltrasound machine witli a software-programmable subsystem. 

15 Tissue vibration imaging can provide additional diagnostic information that is 

currently not available to the clinician using conventional tools. An ultrasound device 
with tissue vibration imaging capability can become a usefiil screening and diagnostic 
tool for paramedics and trauma centers for rapid diagnosis and localization of active 
internal bleeding. Such real-time localization of bleeding can tlien be employed for 

20 targeting and evaluating tlie effects of HIFU and otlier therapies intended to stop the 
bleeding. 

Although the present inveixtion has been described in connection with the 
prefen-ed form of practicing it and modifications thereto, those of ordinary skill in the art 
v^U vmderstand that many other modifications can be made to the present invention 
25 within the scope of tlie claims tliat follow. Accordingly, it is not intended that the scope 
of the invention in any way be limited by tlie above description, but instead be 
deterrriined entirely by reference to the claims that follow. 



